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We study the orthogonality catastrophe due to a parametric change of the single-particle (mean 
field) Hamiltonian of an ergodic system. The Hamiltonian is modeled by a suitable random matrix 
ensemble. We show that the overlap between the original and the parametrically modified many- 

2 

body ground states, S, taken as Slater determinants, decreases like n , where n is the number 
of electrons in the systems, k is a numerical constant of the order of one, and x is the deforma- 
tion measured in units of the typical distance between anticrossings. We show that the statistical 
fluctuations of 5* are largely due to properties of the levels near the Fermi energy. 

PACS numbers: 73.23.-b, 71.10.-w, 73.23.Hk 



I. INTRODUCTION 

The Anderson orthogonality catastrophe (AOC), in- 
troduced by Anderson in 1967tl, is a fundamental ef- 
fect observed in many-body systems. The original work 
addressed the ground state of a finite system consist- 
ing of N non-interacting electrons. Upon the introduc- 
tion of a localized perturbation, this ground state gets 
modified. Anderson has shown that the overlap between 
the original and the modified A^-electron ground states, 
(^jvI^jv)' i s proportional to a negative power of N, and 
vanishes in the thermodynamic limit, hence the catas- 
trophe. Variants of the AOC are at the basis of some 
central themes in solid state physics, including the x-ray 
edge singularity, zero-bias anomalies in disordered sys- 
tems and tunneling into quantum Hall systems. 

The applicability of this concept and attempts to ex- 
tend it to more generic circumstances have been at the 
focus of attention for more than three decades. Partic- 
ularly appealing is the application of AOC ideas to the 
field of mesoscopics. The study of mesoscopic systems 
involves finite size systems where it is usually important 
to account for the dynamics and the thermodynamics of 
the electrons on a quantum-mechanical level. Important 
ingredients in characterizing mesoscopic electronic sys- 
tems include the strength of the ambient disorder poten- 
tial, the system's size and shape and the strength of the 
electron-electron interaction. Evidently, in most cases 
such systems are too complex to be described or ana- 
lyzed exactly. In such situations one needs to resort to 
various approximation schemes. 

Finite size ("zero-dimensional") conductors, a.k.a. 
quantum dots (QDs) have received special attention in 
recent years, partly due to the rich physics-involved, but 
also due to their experimental accessibilityacl. The sim- 
plest scheme to account for a QD with interacting elec- 



trons is the constant interaction modelO. The latter im- 
plies that the system Hamiltonian is given by 

H = H + Hci, (1) 

where Hq is the single-particle Hamiltonian and the in- 
teraction, represented by an infinite wavelength (zero 
mode), is given by 

He i = ^(N - iVo) 2 . (2) 

Here C is the total effective capacitance of the dot, N is 
the particle number operator, and Nq represents a tun- 
able background charge (related to the gate voltage) . For 
a complex QD with diffusive disorder or with chaotic dy- 
namics, the single-particle energy spectrum and eigen- 
functions of Ho should be the subject of a statistical de- 
scription. It turns out that for such systems spectral 
correlations within energy windows up to the Thouless 
energy, -&rh, are well described by the random matrix 
theory (RMT), see for instance Refs. |]j7| The individual 
wavefunctions, hence spectral correlations, remain un- 
changed as we add/remove electrons from the system. 
Other than the Coulomb gap, we do not expect any sig- 
nature of the AOC. 

For various reasons the constant interaction model 
does not provide a satisfactory framework to account 
for key phenomena. The latter include some features 
of the addition spectrum, Coulomb peak-height correla- 
tions and electron scrambling (see e.g. Ref. ||). To im- 
prove on that model one can employ the best effective 
single-particle approach, the Hartree-Fock (HF) approx- 
imation. Expressed in the basis of the exact eigenstates 
of Hq, denoted by {i/> a } with the corresponding eigenen- 
ergies {e a }, the HF Hamiltonian reads 
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-ffHF = ^ S a it 



a, S3 



(3) 



Here v a p a p are the antisymmetrized matrix elements of 
the interaction and h a is the number operator of the 
state i/> a . Note that -Hhf includes only diagonal inter- 
action matrix elements (in the exact eigenstate basis); 
for a short-ranged interaction it has been shown that-pff-. 
diagonal matrix elements are parametrically smalLETtEl 
It is clear that within this approximation the effective 
single-particle states (but not the HF energies!) are in- 
dependent of the occupations of these states, and are 
unchanged upon the addition/removal of electrons hops, 
the system. This is in line with the Koopmans' picture.til 
Only when the single-particle states of the systern-are 
modified upon the introduction of a "perturbation" t3 is 
it meaningful to address the question of the AOC. Here 
the perturbation can be understood as varying a localized 
gate voltageEJ, introducing a static impurity or adding an 
electron to the (N + 1)— st state. The latter motivates 
the study of 



(4) 



where \^n) is the TV-electron system ground-state and 
cjv+i is the corresponding electron creation operator at 
its first single-particle excited state. In all these examples 
(including the HF scheme) both ground states are given 
by Slater determinants. Indeed the Anderson formalism 
deals with the overlap of two Slater determinants. More 
specifically, one of the applications we have in mind is 
the role of AOC in the electronic transport through com- 
plex quantum dots. At very low temperatures their con- 
ductance peak heights in the Coulomb blockade regime 
can be evaluated from the tunneling rates to and from 
the corresjjictfiding many-body ground state of the elec- 
tron islandtil. In the simplest approximation, these rates 
are given by the overlaps between a single-particle wave 
function in the dot and the channel wave functions. How- 
ever, if interactions are taken into account, even in a 
HF mean field approximation, a many-body contribu- 
tion to the tunneling rates has to be considered. A pri- 
ori the inclusion of an additional electron to the island 
may change each individual single-particle orbital negli- 
gibly; however, if the number of electrons in the dot is 
large enough, the new many-body ground state can be 
almost orthogonal to the old one. Some aspects of the 
AOC in the presence of disorder have been previously 
considered, see e.g. Ref. ^|. In the present paper we try 
to circumvent the task of analyzing the AOC in the pres- 
ence of disorder potential, and replace this challenge be 
studying the AOC within the framework RMT. Here the 
role of an added perturbation is played by a parameter 
which is varied. We end up studying ±he AOC in the 
context of parametric random matricesEj. As has been 
argued recentlytLZl this is a good model to describe elec- 
tron scrambling in quantum dots embedding the discrete 
electron number in the dot in a continuous variation of a 
parameter .113 



The paper is organized as follows. In Section || we 
rederive the expression for the ground state overlaps, re- 
viewing the main features of the theory and introduc- 
ing the notation employed throughout this work. Sec- 
tion III is devoted to the presentation of the random 
matrix model used in this study. Section [V contains 



the main body of our analysis, presenting our analytical 
and numerical results for the averaged ground state over- 
laps. The distribution of \S\ is the subject of Section [v|. 
Section VI contains some brief concluding remarks. We 
have also included three appendices. Appendix A quan- 
tifies the accuracy of the "unit volume approximation" 
used by Anderson in his original worklll. In Appendix B 
we discuss some implications of evaluating the ensemble- 
averaged value of the Anderson integral subject to grand 
canonical constraints. Finally in Appendix C we discuss 
the precision of the first order perturbation theory used 
to evaluate Anderson's integral, defined in Section O. 



II. REMINDER ON THE ANDERSON 
ORTHOGONALITY CATASTROPHE: A 
PARAMETRIC APPROACH 



Let us consider systems of electrons whose Hamilto- 
nian, H(X) is a function of a parameter X. In this 
study the generic parameter X can be either continuous 
or discrete. In the case of modelling the change in an ex- 
ternal magnetic field or in the electrostatic potential, X 
is continuous. One can also imagine X as modelling the 
change in the effective single-particle potential due to the 
sequential addition of electrons to the system, in which 
case X is discrete (in the latter case it is necessary first 
to subtract all systematic changes of Jthe single-particle 
spectrum.) It has been shown recentlyIHa that scrambling 
due to the sequential addition of electrons to the QD can 
be embedded in a parametric RM process only if interac- 
tion matrix elements due to the accumulation of surface 
charge are taken into account. The subject of the present 
paper is the overlap S between ground states of systems 
corresponding to H(X) and H (X + 5X). In particular, 
this section is devoted to the presentation of different 
levels of approximation for log \ S\. 

Since we are dealing with non- interacting electrons, 
the ground state wave functions are Slater determinants. 
The overlap S can then be expressed in terms of the 
occupied single-particle wave functions overlap. Let us 
denote by tpk(X) and ipk{X + SX) the (single-particle) 
eigenstates of H(X) and H(X + SX) respectively. We 
define the unitary overlap matrix A as 



A; 



{MXM 3 (x + sx)) 



(5) 



Thus the overlap of the ground states S can be written 
as 



S = dct A° 



(6) 



2 



where the superscript "o" stands for occupied states. Ac- 
cordingly A°° corresponds to the subspace of occupied 
states of the matrix A defined in Eq. (^). 

The catastrophe is manifest by the fast suppression of 
the overlap S as the number of fermions n in the system, 
is increased for a fixed perturbation strength. AndersonEJ 
singled out two basic reasons that make the absolute 
value of the overlap S smaller than one: The rows of 
A°° have norms smaller than unit and they do not form 
an orthogonal set. It is convenient to separate these two 
contributions by introducing a new matrix A°° with nor- 
malized rows, defined as 



=(A°°) ij M 



where the normalization factor Afi reads 



Afi = 



\ 



E 



\S\ is then rewritten as 

n 

\S\ = \detA^\ x 



(7) 



(8) 



(9) 



The absolute value of det A°° can be geometrically inter- 
preted as the volume of a n-dimensional parallelepiped 
defined by a set of non-orthogonal vectors of unit length. 
In Ref . |l| it is claimed that | det A°° | does not depend on 
n and it is close to one. In Appendix [A] we show pertur- 
batively that indeed the most important contributions to 
15*1 come from J^A/i, the deviations from the "unit vol- 
ume approximation" being of higher order in 5X. Aiming 
at obtaining an upper bound of the overlap S we shall 
neglect below the contribution of | det to \S\. We 
thus restrict our analysis to the product of normalization 
factors. We can therefore write 



\S\<Y[Af i = l[(l-P i ) 



1/2 



i=l i=l 



where we have introduced 



Pi= \(MxM(x + sx))\' 



(10) 



(11) 



j=n+l 



Pi measures the probability of the state i, whose energy 
Si is smaller then the Fermi energy, to be spread over 
the ^-components lying above the Fermi surface once the 
Hamiltonian is changed. The product of the normaliza- 
tion factors, Eq. ( |l0| ) , is the first level of approximation to 
the overlap determinant considered here. For the sake of 
convenience, hereafter we shall deal with log | | instead 
of 151 itself, and define 



Eq. (|lj) is still not suitable for an insightful analysis. 
For this purpose a higher level of approximation for S 
was introduced!!!. For small SX the "probabilities" Pi 
will also be small. Thus, by expanding J norm in a power 



series 



1 " P 2 

'norm — ^ / A * 2 
i=l 



(13) 



and retaining only its first expansion term one writes 



1 - 



(14) 



By explicitly expressing / in terms of the single-particle 
overlaps, one arrives at 



n oo 



(15) 



z—l j—n+1 



This is the well known Anderson integral term. The ap- 
proximations /norm and I for log IS] will be analyzed at 



length in Section IV where we present numerical and an- 



alytical results for the random matrix model introduced 
in the forthcoming section. 



III. THE RANDOM MATRIX MODEL 

In this Section we present a model for the statistical 
study of the Anderson orthogonality catastrophe in bal- 
listic ergodic systems. The setting is the same as in the 
foregoing Section. We model the single-particle Hamil- 
tonian H which depends on a parameter X by a suit- 
able assemble of random matrices. In one of its simplest 
formala, such model is realized by 



H(X) = H cos X + U sin X 



(16) 



where Hq and U are two independent matrices of dimen- 
sion N x N, both members of the same Gaussian Ensem- 
ble. In this study, we restrict our analysis to the Gaus- 
sian Orthogonal Ensemble Q3 = 1) and to the Gaussian 
Unitary Ensemble ((3 = 2). Whereas the orthogonal en- 
semble is used to model systems presenting time-reversal 
symmetry, the unitary ensemble models systems where 
the latter sjawnetry is absent. Thus the distribution of 
H and U is 



P( J ff ,C/) = C /3Jv exp 



2N 



(trifoi^+trt/T/t) 



(17) 



-£>g(l-i>0> log \S\ 



(22) where CpN is a normalization constant. As a conse- 
quence of this choice, the matrix elements of H{X) are 
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also Gaussian distributed with zero mean and variance 
parameterized as 

= (SwStf + Sfn&wSiv) ^ (18) 

independent of X. Here the symbols (• • •) stand for en- 
semble averaging. 

For any arbitrary value of X, the resulting mean level 
density Js given by the well-known Wigner semicircle 
formulaB 



(19) 



where e is the single-particle energy. Evidently, the mean 
level spacing at the center of the spectrum (e = 0) is 
A = ttX/N. Accordingly, the average band width is 
E max — £min = 4A (this is in line with standard solid state 
picture whereby the band width is fixed, independent of 
the system's size, while the mean level spacing scales as 
TV- 1 ). 

The typical scale for the parametric variation X* rep- 
resenting the average distance between anticrossings, can 
in generaLb^jCharacterized by the level velocity standard 
deviatiorMy 



A 

X* 



\ 



dX 



(20) 



where £ M (X) is an energy level of H(X) close to the cen- 
ter of the band. In the model defined by Eqs. (16) and 
(Il7|) the average level velocity is zero and 



X* 



A 



V27PV)A 




(21) 



Now it is possible to quantify the effect of 8X in the 
single-particle spectrum by the scaled parameter x = 

sx/x*. 

It remains to specify the many-body part of the model. 
The ground state configuration is generated by popu- 
lating the single-particle levels with n = N/2 fcrmions 
(the spin degree of freedom is not considered here). For 
a given X, the ground state is the Slater determinant 
made up of the lowest N/2 eigenstates of H(X). Now 
the single particle Hilbert space is finite and of size N, 
and the number of particles is N/2. The Fermi level, 
£f = (£n/2+i + £at/2)/2, lies on average at ef — 0. It is 
important to stress that the sums in Eqs. (|l|) and 
have to be modified accordingly. 

Numerical simulations of the model presented above 
are straightforward to implement. We generate a pair of 
matrices H and U of dimension N whose elements are 
Gaussian distributed with zero mean and variance given 
by Eq. (|l8|). Fixing X and x we calculate the eigenval- 
ues and eigenfunctions of the corresponding H(X) and 
H(X + SX). Having computed the complete set of eigen- 
functions, we calculate the overlap matrix A°° and its 



determinant, Sn(X, x). In order to keep the notation 
simple, we shall write down the arguments of S explic- 
itly only when necessary. 

As a guide for the discussions to come let us consider 
a representative member of the ensemble defined by Eq. 
([[!) for ft = 1 and N — 50. In Fig. we show the 
parametric spectrum as a function of X around e = 
(top panel) and the corresponding ground state overlaps 
|5(X)| for x = 0.4 (bottom panel), both obtained by a 
direct numerical procedure. The most striking feature of 
this plot is that large fluctuations in 1*91 are closely cor- 
related with the occurrence of narrow avoided crossings 
at the Fermi surface. In other words, to each narrow gap 
(avoided crossing) there corresponds a small overlap \S\. 
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FIG. 1. Top panel: typical energy levels as a function of X. 
The filled region is the gap between the last occupied and the 
first empty single-particle levels. Bottom: the three curves 
represent the exact overlap S, exp(/ n orm), and exp(J) (see 
text for definitions) as a function of X. They respect the or- 
dering 5* < exp(7 n0 rm) < exp(J). The parametrical distance 
is set to x = 0.4 (5X ~ 0.13) and the dimension of the single 
particle Hilbert space is N = 50. 

The figure also suggests that the approximation schemes 
presented in Section |l| fail in the vicinity of narrow gaps. 
These observations indicate that an accurate analytical 
estimate of the average \S\ and its variance requires a 
good handle of the "two- level problem" of narrow avoided 
crossings. Had we considered integrable (or mixed) sys- 
tems, this issue would become even more important due 
to the absence (or suppression) of level repulsion. 
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IV. AVERAGE GROUND STATES OVERLAP 

As we shall argue in the next section, the distribution 
of S can be anomalously broad. As a first step though, we 
shall consider in this section ensemble averaged quanti- 
ties. When the distribution is narrow, as is often the case 
with thermodynamic quantities, one may interchange 



The occurrence of small gaps at the Fermi surface, 



(log|S|)«log<|S|>. 



(22) 



Evidently, when in doubt, the meaningful quantity is the 
r.h.s. of the equation above. Technically, though, the 
more accessible quantity is the l.h.s. of Eq. (22). This 



section is devoted to the statistical study of upper bounds 
for (log LSI) within the approximation levels presented in 
Sect ion O and the random matrix model described in Sec- 
tion [II]] Our analytical results are complemented with 
numerical simulations. 

In Eq. Jl5| ) the ground states overlap can be easily 
evaluated by considering the single-particle overlaps in 
first order perturbation theory, namely 



5 Hi, 



(23) 



Here 6H w (— HqsIxiX + U cosX) SX and {ej} are the 
eigenvalues of H(X). The matrix 5H and the set of 
eigenvalues {sj} are statistically independent due to the 
invariance of the considered ensembles under orthogonal 
(J3 = 1) or unitary (fl = 2) transformations. Recalling 
Eq. ( aa ) (log 1 5 1) now reads 



N/2 n 



(logLS|><<7> = -i]r £ 



1=1 j=N/2+l x v J ' 



(24) 



In Appendix ^ it is shown that the corrections responsi- 
ble for the inequality in Eq. (p4| ) are of fourth order in a 
perturbation expansion. Note that only off-diagonal SHy 
matrix elements contribute in Eq. (g3() since the i-states 
lie below the Fermi surface and the label j corresponds 
to states above it. The ensemble average over \5Hij\ 2 , 
defined in Eq. (pit), yields 



N/2 N 
t=l j=N/2+l 



A 2 



(25) 



where we have used X* of Eq. ( |2l| ) to express ( log |5|) 
in terms of the dimensionless variable x = SX/X*. 

A satisfactory accuracy of the approximation intro- 
duced by Eq. ( p3| ) requires SHij/(ej — £j) to be small. 
Such a condition is translated to (I) by examining sepa- 
rately x 2 and (A 2 /(sj — Si) 2 ). Viewed from the perspec- 
tive of the parametric framework, the approximation is 
under control since the physical situations in_mind, see 
Section ffl, call for x <C 1, or at least for x < 10. By con- 



trast, the average over A /(sj 
discussion. 



Si) requires a careful 



ej - Si < A for j 



1, causes the breakdown of the 



perturbation approach employed in Eq. H23J) , to calculate 
(tpi(X)\ipi+i(X + SX)). In that case even a small varia- 
tion of x <C 1 may give rise to a significant mixing of an 
originally occupied level (the last one) and an empty one 
(the lowest original vacant level). Such a situation calls 
for a non-perturbative solution. Indeed it is tempting to 
justify the approximation used in Eq. (|2|) invoking the 
presence of level repulsion: the large errors introduced in 
the calculation of \S\ at small gaps are minimized by the 
rareness of such events. Unfortunately, such a scheme is 
doomed to fail due to the following reason: For neigh- 
boring levels 



A 2 



i+1 



dss- 2 Pp(s) , 



(26) 



where Pp(s) is the nearest neighbor spacing distributioni. 
Since, Pp(s) ~ s 13 for s <g; 1 the average on the l.h.s. of 
Eq. (^6|) diverges for the orthogonal symmetry (although 
not for the unitary case). This motivates the employ- 
ment of a non-perturbative approach to account for the 
effect of the "crust levels" (near the Fermi energy) on the 
ensemble-averaged Anderson integral. With this proviso 
in mind it should be also realized that j — i = 1 corre- 
sponds to a single term in the double sum of Eq. ( p5|) . As 
concluded from Fig. |l| this term has a major contribution 
to the large fluctuations in \S\ and will be the subject of 
analysis of Section [y| By now it is only necessary to an- 
ticipate that the occurrence of narrow gaps contributes 
to (log IS* |) with an additional factor which depends on 
x and on the size of the narrow gap. Note that the pro- 
cedure presented in Eq. ( p6| ) represents averaging under 
canonical conditions, namely averaging over systems with 
a given particle number AE3 or averaging over both im- 
purity realizations and N. Such procedures are referred 
to as strong and weak canonical averaging, cf. Ref. ||^. 
There are experimental setups where it is the chemical 
potential, fi, which is the controlled parameter. In such 
circumstances it is more appropriate to employ a grand- 
canonical averaging procedure. This is briefly discussed 
in Appendix |^. 

As the value of j — i becomes larger the perturbative 
approach works increasingly better. Moreover, at the 
same time the fluctuations of Ej — Si relative to (j — i) A 
decrease as a consequence of the spectral rigidity. These 
matters are discussed in Appendix ^|. In this regime we 
do not expect to introduce a large error (independent of 
N) replacing the spectrum of H(X) by its average spec- 
trum, that is 



(27) 



(Ej-eif/ ({ej)-{ei)f 



The latter approximation allows us write 
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-A/2 



2A 



2A 



A/2 



P{£j)p{£j) 



(28) 



with the mean level density p(e) given by the Wigner 
semicircle law, Eq. (|l|). Changing variables and defin- 
ing a = 7r/(4JV) we arrive at 



(I) =~ 1 x 



2 1 du I dv - ; . (29) 



a 



(u + v) 2 



We can now isolate the singularity at a = (N large) 
and rewrite the integral on the r.h.s. of Eq. (p9| ) 



i ,i 
du / aw 



1 



i 



(u + v) 2 



(u + vf 



(30) 



While it is the first term under the integral in Eq. ( |30| ) 
which is ultimately responsible for the catastrophe, the 
second one gives only a constant as a — > 0. The final 
result is therefore 



(I)=-L 2 (\ogN + C) , 



(31) 



where C — — log(7r/2) — 7r 2 /8. The latter constant has to 
be taken with caution due to the approximations made. 

Summarizing our results, we have several levels of ap- 
proximation to the exact value of (log|5|). The product 
of the normalization factors, which neglects corrections 
to the volume of the parallelepiped, assumed to be unity, 
yields 



N/2 



(log|S|)«(7 norm ) = i^(log(l-P i )) 



(32) 



By keeping only the first term of the expansion of the 
logarithm in I n orm we have 



N/2 



(\og\S\) ^ (I) ^ -^J2(P t 



(33) 



By calculating I using first order perturbation theory and 
a smoothed spectrum, we obtain the analytical estimate 
(7), given by Eq. ©. 

For small values of x these quantities are close to each 
other. They are ordered as follows: 



(log |^|) < (7 norm ) < (I) < (I) 



(34) 



We turn now to our numerical analysis. For different 
values of x and N we have evaluated (log|S|), (/„), (I), 
and (I) ensemble averaging over M realizations of H(X) 
which is defined in Section III. Each simulation is per- 
formed at the cost of the order of Mx (fix N) 3 operations, 
imposing a computational constraint on the procedure. 
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FIG. 2. The average log|5(x)| as a function of TV. The 
number of realizations for each N is M = 10 4 and (a) x = 0.1, 
(b) x = 0.2, (c) x — 0.5, and (d) x — 1.0. Open dots stand 
for the exact (log \S\), filled dots for (/norm), squares for (I), 
whereas the solid lines represent (7) . 

In Fig. U we show a comparison between the differ- 
ent approximations for (log|5|) as a function of N for 
(3 = 1. We chose four representative values of x and fixed 
71/ = 10 4 . We observe that in all approximation schemes 
(log \ S\) displays the slope predicted by Eq. (|3l|). In view 
of the large sample sizes, the fluctuations in the mean val- 
ues indicate the corresponding distributions are charac- 
terized by large standard deviations, as we discuss in the 
following section. As analyzed in Appendix ^| the first 
order perturbation theory estimate breaks down when 
used for levels in the vicinity of the Fermi surface and/or 
for not sufficiently small values of i. As the latter is 
increased we even expect deviations from the predicted 
slopes for (I) versus log/V. In Fig. ^ such discrepan- 
cies are only noticeable for x = 1 at N < 50. We con- 
clude that the power law suppression of |5| as a function 
of N is very robust. The discussion presented in Ap- 
pendix |c| suggests why our theoretical estimate does not 
always serve as an upper bound for all our approxima- 
tion schemes. This is done by noting two facts, (i) For 
small values of x the perturbation approach underesti- 
mates (|^4jv/2,A72+j| 2 )j hence overestimates the overlap. 
In such cases we are guaranteed an upper bound, (ii) 
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For larger values of x the situation is quite the opposite. 
When this happens the constant C does not provide any 
longer an upper bound. Results for the (3 = 2 case behave 
in the very same way as for (3 = 1, with slopes following 
Eq. (31). (We have thus decided to omit a corresponding 
figure for the unitary case). 

Figure || also confirms Anderson's claim that the "unit 
volume" corrections to (log 151) are small, at least for 
x <C 1, as indeed is seen in panels (a) and (b). While 
within our statistical precision we do not observe any N 
dependence for the unit volume corrections, the latter 
cannot be ruled out. 



V. DISTRIBUTION OF OVERLAPS AND THEIR 
LARGE FLUCTUATIONS 



This section is devoted to the analysis of the distribu- 
tion of the overlaps between ground states, P(|5|). We 
show that the large fluctuations of S(X) occurring in the 
vicinity of small gaps at the Fermi surface are well de- 
scribed by a 2 x 2 model, otherwise exp((/)) provides a 
satisfactory estimate for \S\. 

Let us start quantifying the influence of the narrow 
avoided crossings on the overlap determinant. In the 
proximity of a narrow gap the single-particle overlap 
(ipN/2(X)\ipN/2(X + 5X)) will be small and rapidly vary- 
ing with X. Consequently, as it was qualitatively estab- 
lished in Section ||, it will dominate the fluctuations of 
the single-particle overlap matrix determinant. In this 
situation we expect that the behavior of S will be cap- 
tured by the approximation 

| S N (X, SX)\^K N (X, SX) | {^ N/2 (X) | i, N/2 (X + 6X))\ . 

(35) 

This approximation is just the first term of the expan- 
sion of the determinant of A°° along its iV/2-th row. In 
other words, we calculate |5| as if all fluctuations are due 
to the interaction of the highest occupied single-particle 
level with the lowest empty one. The factor K accounts 
for the contribution of all remaining states and will be 
approximately constant: 



K N (X,6X) ^ (\S N ^{X,SX)\) 



(36) 



The results presented in Fig. || confirm that this approx- 
imation works impressively well in the vicinity of narrow 
gaps and for small values of x. The system is the same 
as that in Fig. [l], but smaller deformation is considered, 
x = 0.2, as well as a smaller parameter range for X. 
The latter contains, in the present example, two narrow 
gaps (cf . Fig. [j]) . The small shift between the two curves 
can be attributed to the K term. The proposed approxi- 
mation gives rise to spurious peaks whenever an avoided 
crossing (narrow gap) between the occupied levels N/2 
and N/2 — 1 is encountered. Such cases are beyond the 



scope of our approximation: the factorization put for- 
ward by Eq. (B5 



is no longer valid, and the gap at the 
Fermi energy is certainly not small. The spurious peak 
atl« 1.45 in Fig. pi represents such a situation, as can 
be verified from Fig.lU 



x 

00 



oq 
d 



CO 

d 



d 




FIG. 3. Overlaps of Slater determinants for x = 0.2 as a 
function of X for N = 50. The exact determinant is shown as 
a black line. The grey line corresponds to the approximation 
\S\ « \(i> N/2 (X)\ip N/2 (X + 5X))\. 

We have thus demonstrated heuristically that small 
Anderson overlaps arise due avoided crossings at the 
Fermi level, involving the highest occupied and the lowest 
unoccupied levels. This observation leads us to conjec- 
ture that 



P(\S\) & P{K ■ \S\) for \S\ < 1 



(37) 



where S = (^ N / 2 {X)\^ N / 2 {X + 5X)). The deriva- 
tion of an expression for P(|5|), which calls for a non- 
perturbative approach, is discussed now. Due to level 
repulsion two levels rarely come very close to each other. 
When this happens, it is possible to effectively model a 
narrow avoided crossing by a 2 x 2 matrix, since it is quite 
unlikely to have yet another level very close by. This 
model also rids of the undesired occurrence of spurious 
peaks as the one shown by Fig. 0: they cancel between 
the factor K and S. While this is quite a simple model, 
the calculation of the distribution -P(|<SU involves some 
cumbersome multidimensional integralscj. 

For the orthogonal case, (3 = 1, the integrations can be 
simplified through a geometrical construct, which makes 
it possible to obtain P(|5|) analytically within a 2 x 2 
parametric random matrix model. Let us write our model 
Hamiltonian, Eq. (pif), H(X) = if cosA + UsinX as 



H(X) = cosX(H a + Ut&nX) 



(38) 



Let us also parameterize its eigenvector as ip 
(cos t9, sin (9). Then, by defining the vector h(X) 
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(2_ffi2, H22 — Hix), it is straightforward to show that the 
eigenvector equation can be written as 



(2H 12 ,H 22 - H n ) • (cos 29, sin 29) = , 



(39) 



so that the eigenvectors are determined solely by the vec- 
tor h(X). Moreover, the angle between ip(X + SX) and 
ip(X) is half the angle a between h(X + SX) and h(X), 
that is |5| = cos a/2. Now the problem is reduced to 
finding the distribution of a. Let us set X = and in- 
troduce 

h(X = 0) = h = (2[tf ]i2, [H } 2 2 - {H ]n) 

h(SX) = h t = (2[flo]l2, [H Q ]22 - [Bo]ii) + 

{2U 12 ,U 22 -U n )t, (40) 

with t = tan<W. The usefulness of this geometrical 
construction becomes clear now. We use the fact that 
h t = ho + ut and 



hf h = \h t \\h \ cos a , 



(41) 



where the vectors h t and ho can be expressed in terms 
of \ho\, I it I, 9 U , and a. The integration of P(ho,h t ) over 
ho, u and 9 U , readily gives P (a). 

The distribution P(ho,h t ) is written in terms of Hq 
and U matrix elements, Eq. (O), as 



PUin.h,) = ! 2 J''{ h ,U=j[h t 



hn 



(42) 



The integration over the vectors ho and h t , keeping their 
relative orientation fixed, yields 



P(a) 



2t 2 

— / 1 

I" JO 



(1 + t 2 COS 2 



2 cos 4> sin cos a) 2 
(43) 

The latter is expressed in a closed form as 
t 2 



V ' Tr(sm 2 a + t 2 ) 



[1 + G(a)] , 



(44) 



with 

G(a) 



\J sin 2 a + t 2 



7T 

— + tan 



cos a 



vsh 



t 2 



(45) 



Since we are not really interested in P(a), but rather in 
P (\S\) instead, a last change of variables is needed to 
arrive at the main result of this section 



P(\S\) 



P[a = 2 cos" 



1-\S\< 



(46) 



In particular, the probability density of having a null 
overlap is non- vanishing and is given by 



P(\S\=0) = 



1 - 



SX 
tan SX 



(47) 



Comparison between the analytic form of P(K ■ \S\), 
Eq. ( [i"6"| ) [cf. also Eq. (p7|)], and our numerical study of 
the exact distribution of overlaps -Pd^l) is displayed in 
Fig. ^. It is evident that the 2-level picture reproduces 
the tails of the distributions P(|5|) for large N remark- 
ably well, provided we use the same relative deformation 
x. Note that for N = 2 and (3 = 1, one has X* = ^/tt, 
differing from the large N limit given by Eq. (|2l"|), see 
for instance Ref. ^5|. The correction due to K^, a rescal- 
ing of the horizontal axis, is negligible in this case (from 
Fig. H one finds that K ps 0.98) and was not included. 




0,0 0,2 0,4 0,6 

ISI 



1,0 0,0 0,2 0,4 0,6 0,8 1,0 
ISI 



FIG. 4. Log-normal plot of the distribution of overlaps 
Pd^l). We compare the analytical prediction based on the 
case N = 2 with numerical simulations over 10 pairs (Ho, U) 
(histograms), for TV = 100. The relative deformation is (a) 
x = 0.1 and (b) x = 0.2. 



The agreement between the output of our 2x2 model 
and the exact diagonalization becomes even more evi- 
dent by analyzing the cumulative overlap distribution in 
a log-log plot, as is seen from Fig. || 
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100 



estimatedBHa to be 



1E-4, 



FIG. 5. Log-log plot of the cumulative distribution of over- 
laps N(\S\) = / |S| d\S\ P(\S\). We compare (i) the analytical 
prediction based on the case N = 2 and (ii) numerical sim- 
ulations over 10 pairs (Hq,U) (dots), for N = 50. In both 
cases the relative deformation is x — 0.2. 



VI. CONCLUDING REMARKS 

We have studied the orthogonality catastrophe due to 
a parametric change of the single-particle of an ergodic 
system. The Hamiltonian was modelled by a suitable ran- 
dom matrix ensemble. We show that the average over- 
lap between the original and the parametrically modified 
many-body ground states, taken as Slater determinants, 
decreases like n~@ x / 4 , where n is the number of electrons 
in the system and x is the deformation measured in units 
of the typical distance between anticrossings. We have 
also shown that the fluctuations of log l^l are enormous. 
To account for the latter in the orthogonal case (f3 = 1), 
we have put forward a simple 2x2 matrix model and 
employed it to obtain a prediction for -Pd-SI) for |S| <C 1, 
in good agreement with our numerical analysis. In the 
unitary case (/3 = 2) the fluctuations are smaller, but still 
quite significant. Here also it was showncilthat the 2x2 
model works well, but no simple analytical expression is 
avaliable. 

This model study constitutes a first step towards un- 
derstanding the relevance of the orthogonality catastro- 
phe for ballistic ergodic systems. One improvement to 



the theory presented in Section IV should arise from the 
fact that for a generic dynamical system the "pertur- 
bation" matrix, SHij , will in general have a finite band- 
width b. How are our results changed in this case? There 
is extensive literature dealing with the matrix elements 
properties of low-dimensional dynamical systems. In par- 
ticular, for the case of a two-dimensional chaotic ballistic 
system, the bandwidth for a "generic" perturbation was 



tA 



(48) 



where r is a time needed for a particle to traverse the 
system, i.e., 

L 



(49) 



Here L is the linear dimension of the system and vf is the 
Fermi velocity. Relating the 2D mean level spacing to the 
system's size, A = 27r?i 2 /(mL 2 ), and using ef ~ nA, we 
obtain a simple relation between the bandwidth and the 
number of electrons, namely, b ~ ^fn. In line with the 
present study, we assume that within the band — < b, 
the fluctuations of SH.^ are statistically well described 
by RMT. Accordingly, the average of log \ S\, obtained in 
Eqs. @ , (H), becomes 



(log|5|) « -|x 2 log. 



(50) 



We stress that in order to write Eq. (|50|), two impor- 
tant assumptions have been-jpade. (i) The Hamiltonian 
H(X) must be fully chaoticEH (ii) The perturbation must 
be generic in the sense defined in Ref. ^8|. If such con- 
ditions are not met, the statistical distribution of SHij 
could strongly depend on \j — i\, for \j — i\ < b. 

The study of the overlap between of many-body states 
subjected to a perturbation has also been a traditional 
subject of investigation in nuclear physics problemscj. 
More recently there is a renewed interest in examining 
overlaps among excited statesS3. This study is much 
more involved than ours, since the states considered are 
superpositions of Slater determinants, rather than a sin- 
gle one as studied here. A full theory for such a case is 
still lacking; available numerical evidence (evidently for 
small systems) indicates that the scaling of the overlaps 
with n follows a Gaussian, a fact which is yet to be ex- 
plained. 
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APPENDIX A: A PERTURBATION SERIES FOR 

logS 

In Section |n| we showed that the overlap of the ground 
states S can be written as S — detA 00 , with the matrix 







elements of A defined by Eq. (||). This is done for the 
case where time- reversal symmetry is present (absent), 
namely the orthogonal (unitary) symmetry. The nota- 
tion "oo" stands for "occupied-occupied" states (in this 
Appendix we shall also use "oe" for "occupied-empty" ) . 
Within the Rayleigh-Schrodinger perturbation theory the 
overlap matrix is expanded as 



A = 1 + eB + e 2 C + e 3 D 



(Al) 



Since AA* = 1, i.e., using the fact that this series is 
normalized, it is implied that 

B + B f = 
C + C ] + BB^ = 
D + D ] + CB ] + BC ] = , (A2) 

etc.. In particular, one knows from perturbation theory 
that the first order correction B is an antihermitian ma- 
trix satisfying Ba — and 



eBij = 1— 

Ej — Ei 



for 



J > i 



(A3) 



To obtain a perturbative series for log |5| we use Eq. ( Al) 
and the identity log det = tr log to write 

log det A°° = tr log (l + eB°° + e 2 C°° + e 3 D 00 + . . .) . 

(A4) 

By expanding the logarithm and regrouping the terms in 
the sum order by order in e, we obtain 

1 9 

-B°° 2 ) + . 



logdetA 00 = tr 



eB°° + e\C° 



The linear term in e vanishes since Bi. 



first two equations in (A2) one arrives at 
1 
2 



(A5) 
0. Using the 



logdetA 00 = ~^e 2 tr (^B oc B oct - C°° + C oot 



0(e 3 ) . 
(A6) 

The contribution from c°° — C° oJ< being purely imaginary 
corresponds to a phase in det A°° , so that 



log | det A°° | = --e 2 trB oc B oct + C(e 3 ) 



(A7) 



If we average over the parametric random matrix ensem- 
ble, Eq. (|l8[), the third order terms disappear: 

1 



(log|detA°° 
or, equivalcntly, 
(log|detA°°|) 



-e z tT(B oc B c 



0(e 4 



\SHi£ 
{Ei-Etf 



(A9) 

where j and i run over the occupied and empty states, 
respectively. 

This calculation also shows that 



exp(log|detA°°|) - (| det A°° |) = C(e 4 ) 



(A10) 



APPENDIX B: GRAND CANONICAL 
AVERAGES 



In this appendix we study a slightly different ques- 
tion, complementary to the ones we have addressed so 
far. It is related to the averaging procedure. Suppose 
we have an ensemble of systems, each of which having 
a given number of electrons, N, which is kept fixed as 
the parameter X is varied. Let us assume, though, that 
the preparation of this ensemble of systems is performed 
grand canonically, i.e., by attaching each sample to a 
weakly coupled particle reservoir at a given chemical po- 
tential /i, equilibrating, and then removing this coupling. 
The number of electrons in each system may vary, but is 
kept fixed during the "measurement" (i.e., varying X). 
This procedure has been defined in Refs. |26|,|23|as grand- 
canonical — canonical. Under these conditions we employ 
a grand-canonical averaging scheme to calculate (log 
Here we expect that large fluctuations due to avoided 
crossings will be suppressed, since the statistical weight 
of an avoided crossing needs to be mOjdified from P(s) 
(canonical) to sP(s) (grand canonical)! 2 -! This reflects 
the fact that the probability to place the chemical po- 
tential in a given gap (between the last occupied and the 
first vacant level) is proportional to the size of the gap. 

The results of our simulations using the grand canon- 
ical ensemble are summarized by Fig. ^. There, as in 
Fig. |^, we show (log 151) and its different approximation 
schemes (explained in Section IV) as a function of N. 



-0,008 
W -0,012 



(A8) S 



-0,016 




x=0.1 



10 



100 



-0.03 



-0,04 



-0,05 



-0,06 



-0,07- 



-0,08 




N 



N 



FIG. 6. The average log|S(x)| as a function of TV. The 
number of realizations for each TV is M = 10 4 and (a) x = 0.1 
and (b) x = 0.2. Open dots stand for the exact (log|5|), 
filled dots for (/norm), squares for (J), whereas the solid lines 
represent (7) . See Section ^ for the definitions. 
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By comparing the results of Fig. g with the ones shown 
above we observe that the values of (log \S\) , as well as 
(/ noral ) and (I), obtained by the canonical simulation are 
systematically reduced with respect to the ones obtained 
using the grand canonical averaging procedure. This is in 
line with the reasoning that the grand canonical averag- 
ing suppresses large fluctuations corresponding to small 
gaps. Thus it eliminates the long small |5| tails of P(|5|) 
characteristic of the canonical ensemble for x <C 1. Wc 
also observe that both ensemble averaging procedures 
tend to render similar results as the value of x is in- 
creased. (For this reason we refrain from showing the 
grand canonical averages corresponding to x = 0.5 and 
x = 1.0 as depicted in Fig. g.) 

Figure [7] quantifies the discussion presented in the fore- 
going paragraph. The object of study is the standard de- 
viation S log \S\ = [((log|S|) 2 )-(log|S|} 2 ] 1/2 - For a fixed 
matrix size N — 50 and M = 10 4 realizations we compute 
the ratio R between the canonical and the grand canoni- 
cal S log \S\ as a function of x. As expected R(x) > 1 for 
all investigated values of x. Furthermore, as we increase 
x the occurrence of small gaps becomes less important 
and the canonical fluctuations get closer to the grand 
canonical ones. 



with the approximation discussed in Section IV 




FIG. 7. The ratio between canonical and grand canonical 
standard deviation 51og|S(a;)| as a function of x for a fixed 
N = 50 and M = 10 4 . 



N/2,N/2+j I 



f3x 2 A 2 



2{ £ N/2+j — SN/2) 



(C2) 



The quantity plotted in Fig. |8| is the ratio R = 

{\A N /2M/2+j\ 2 )/{\A N / 2 ,N/2+ 3 \ 2 )- For j small as com- 
pared to N/2 the average semi-circle spectrum can be 
approximated by a "picket fence" with spacing A, mean- 
ing that (EN/2+3 ~ e N/2) 2 — j 2 A 2 . The ratio R is then 



R{j,x) 



2f 

px 



2 (\^ N/2 (X)\^ N/2+j (X + SX))\ 2 ) 



(C3) 



The figure clearly displays the following features. For x 
fixed, the ratio R —> 1 for large values of j. When x < 1 
the approximation works well except for the few states 
with, say, j < 5. In the case of j = 1, even if x <C 1, per- 
turbation theory breaks down due to very narrow avoided 
crossings. 



ro - 



cm - 



X 




FIG. 8. Test of first order perturbation theory plus av- 
erage spectrum. We plot the ratio R, the quotient "ex- 
act/approximate" (see text), as a function of the energy dif- 
ference j for different values of the parametrical distance 
x = 5X/X* (x = 0.2,0.3,0.4,0.5,0.75,1.0; the curves with 
increasing x correspond to the ones with decreasing values 
for the first abscissa point j = 1). Each curve was calculated 
by averaging over 10000 pairs (Ho, U) with TV = 100. 



APPENDIX C: ACCURACY OF THE 
APPROXIMATIONS FOR THE AVERAGE OF 
SINGLE PARTICLE OVERLAPS 



Here we exhibit a test of the combination of first order 
perturbation theory and an average spectrum to account 
for averages of single particle overlaps. We want to com- 
pare the average 

(\A n/2 ^n/2+j\ 2 ) = (\(iP N /2(XM N/2+j (X + 5X))\ 2 ) 

(CI) 
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